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We extend to additio nal probability me asures and scenarios, certain of the recent results of 



Krattenthaler and Slater (quant-ph/96 12043) — whose original motivation was to obtain quantum 
analogs of seminal work on universal data compression of Clarke and Barron. KS obtained explicit 
formulas for the eigenvalues and eigenvectors of the 2 m x 2 m density matrices derived by averaging 
the m-fold tensor products with themselves of the 2x2 density matrices. The weighting was done 
with respect to a one-parameter (— oo < u < 1) family of probability distributions, all the members 
of which are spherically-symmetric (5'l7(2)-invariant) over the "Bloch sphere" of two- level quantum 
systems. For u = ^, one obtains the normalized volume element of the minimal monotone (Bures) 
metric. In this paper, analyses parallel to those of KS are conducted, based on an alternative 
"natural" measure on the density matrices recently proposed by Zyczkowski, Horodecki, Sanpera, 



and Lewenstein ( quant- ph/9804024). The approaches of KS and that based on ZHSL are found 
to yield [1 + ?J identical 5(7(2) x S'm-Invariant eigenspaces (but not coincident eigenvalues for 
m > 3). Companion results, based on the SU(3) form of the ZHSL measure, are presented for 
the twofold and threefold tensor products of the 3x3 density matrices. In the former case, there 
are six invariant eigenspaces of nine-dimensional Hilbert space, and in the latter, seventeen, of 
twenty-seven dimensional Hilbert space. We find a rather remarkable limiting procedure (selection 
rule) for recovering from these analyses, the (permutationally-symmetrized) multiplets of SU(3) 
constructed from two or three particles. We also analyze the scenarios (all for m = 2) TV = 4, 
TV = 2 x 3, JV = 2x3x2 and TV = 3 x 2 x 2 and, in addition, generalize the ZHSL measure, so that 
it incorporates a family of (symmetric) Dirichlet distributions — rather than simply the uniform 
distribution — defined on the (TV — l)-dimensional simplex of eigenvalues. 

PACS Numbers 03.65.Fd, 03.67.Hk, 02.20.Qs, 05.30.-Ch 
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I. INTRODUCTION 



Wootters, in a paper entitled "Random Quantum States," stated that "there does not seem to be any natural 
measure on the set of all mixed states" . Relatedly, Jones Q asserted that the problem in the Bayesian treatment 
of "the more realistic experimental case of mixed input states . . . lies in selecting a good prior on mixed density 
matrices". Also, somewhat earlier, Band and Park || were concerned with the "lack of a rational axiom of prior 
distribution over the entire domain of density operators" . 

Contrastingly, however, Zyczkowski, Horodecki, Sanpera and Lewenstein [Q have recently proposed "a natural 
measure in the space of density matrices" and used it to "estimate the volume of separable states, providing numerical 
evidence that it decreases exponentially with the dimension of the composite system" . Also, the present author, in a 
series of papers J5|-|| (cf. |l(], sees. 4 and 5] and 0) has pursued the strategy of attempting to normalize the volume 
elements of certain "monotone" Riemannian metrics Jl3| — in particular, the minimal (Bures |lj,|l^]) and maximal 
monotone ones — and use the results for information-theoretic/Bayesian, as well as thermodynamic purposes. In ||, 
it was concluded — supportive of certain earlier results of Petz and Toth fl^| — that among the family of monotone 
metrics, the maximal was most noninformative in character. 

In this communication, we seek to compare and contrast these various measures based on monotone metrics with 
the measures newly presented in Q. In doing this, we immediately encounter a most interesting difference. While 
in the measures are directly defined over the (TV 2 — l)-dimensional convex sets of N x N density matrices, in 

H they are taken to be the products of the (7V 2 -dimensional) Haar (invariant) measure over U(N) and the uniform 
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measure over the (N — l)-dimensional simplex spanned by the N (nonnegative) eigenvalues of the N x N density 
matrices. Thus, in B, the measures are defined in the (larger) (N 2 + N — 1 )-dimensional space. 

II. THE CASE N = 2 
A. The subcase m — 1 

We have been unable to obtain an explicit (degenerate/singular) transformation or mapping to pass from these 
higher-dimensional measures to the lower-dimensional ones, even for the case N = 2. Notwithstanding this, we still 
proceed analytically by, first, averaging the 2x2 density matrices with respect to these differing measures. Such 
matrices are parameterizable, using the conventional "Bloch sphere" representation p™, as, 



1 / 1 + r cos i? r sin d cos $ + ir sin $ sin ( 

2 \r sin •& cos d — ir sin i9 cos $ 1 — r cos i9 



(1) 



(Q < r < 1, 0<?9<7r, < </> < 2tt). Averaging these matrices with respect to any of the (normalized) measures 
used by Slater, expressed in the same variables, we simply obtain, as seems evident, the familiar density matrix of 
the fully mixed state, 



1/2 
1/2 



(2) 



We obtain the same result using the measure of Zyczkowski et al 0, but now expressing ([!]) not in terms of r, 1? and 
(f>, but rather, say, spherical polar coordinates in /owr-dimensions |]l7| , eq. (3.129)]. We, then, perform the averaging 
making use of the associated Haar measure @ cqs. (3.130), (3.131)]. 

B. The subcases m — 2,3 

We enter a less intuitive realm apparently, however, if we average with respect to the same normalized measures 
(that is, probability distributions), the 4x4 density matrices, which are the tensor products of the 2x2 density 
matrices with themselves. (We note that these averaged matrices and all the ones subsequently discussed, describe by 
their very construction, states which are separable, that is, classically correlated p~8[ . For an interesting discussion of 
the significance of the tensor product in quantum computation, see |ll|.) In fact, in (cf. [pl[), Krattenthaler and 
Slater derived a general formula for the entries of 2 m x 2 m density matrices of this type (as well as their eigenvalues 
and eigenvectors), when the averaging was performed with respect to any member of the one-parameter (u) family of 
probability distributions over the Bloch sphere, 



r(| -u)r 2 sintf 
7rir(l - u)(l - r 2 ) 1 



oo < u < 1 (3) 



The asymptotics [m — > oo] of the relative entropy of the 2 m x 2 m m-fold tensor products with respect to their 
corresponding averages, as a function of u, was also obtained in ]2C|], and its relevance to "universal quantum coding" 
discussed (cf. ||J2|]). 

The case u = 4 of @ gives us the (normalized) volume element of the minimal (Bures) monotone metric (cf. |ic| , 
eq. (30)]). The maximal monotone metric corresponds to u = |. Its volume element is not normalizable over the 
entire Bloch sphere. But it is normalizable, if we remove an e-neighborhood of the spherical boundary (r = 1). Then, 
it is possible — proceeding in Cartesian coordinates and taking the limit (e — > 0) of a certain ratio — to obtain two- 
arid one-dimensional (marginal) probability distributions Q . (The one-parameter family ([}]) was employed in [20|] for 
its computational tractability, in addition to its intrinsic interest. This line of work constituted an effort to develop 
quantum analogs of recent seminal results of Clarke and Barron concerning universal data compression [p4| , p5| .) 

Now, let us, for the scenarios m = 2 and 3, average the m-fold tensor products with themselves of the 2x2 density 
matrices, first, with respect to the (normalized) measure utilized by Zyczkowski et al, which is, representing the Eulcr- 
Rodrigues parameters in terms of spherical polar coordinates in four dimensions (0 < 4> < 27r,0 < 6 < ir,0 < \ < 7r), 
0, eqs. (3.129)-(3.131)], proportional to 



3 



sin 2 x sin 9. 



(4) 



The 4x4 and 8x8 results obtained are precisely the same, as those found by averaging the 2 m x 2 m density matrices 
(parameterized as in p)) using the probability distribution for u = — 2 in (Q). For to — 2, this common average is 
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(not simply the diagonal matrix with entries j, as might have been naively anticipated on the basis of (^)) and for 
to = 3, 
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(We observe, somewhat relatedly, that Zyczkowski et al "could not resist the temptation to investigate the mean value 
of [the degree of entanglement] over random density matrices generated" according to their measure [0, App. B].) 



C. Monotone metrics 

So, we are able to conclude that for the case N = 2, the measure employed by Zyczkowski et al does not correspond 
to that for either the minimal or monotone metric. In fact, based on the evidence so far presented, one might 
conjecture that it does not correspond to any monotone metric at all. We say this based on the proposition that for 
any monotone metric one can associate a "Morozova-Chentsov function", expressible as c(x,y) = l/yf(x/y), where 
/(l) = 1 and f(t) = 1 ). For the family of probability distributions (|3|), we have that |l3|, eq. (3.17)] 

f(t) = \±±1L , (7) 

In the minimal (u — i) and maximal (u = |) cases, f(t) = and respectively. Now, both of these functions 
are operator monotone — which constitutes a necessary and sufficient condition for the monotonicity of the associated 
metric. (A function / : R + — > E is called operator monotone if the relation < K < H implies < f(K) < f(H) 
for any matrices K and H of any order. Ordinary monotonicity is obtained if K and H are simply scalar quantities.) 
But for the choice u = — 2, which yields (^) and (||), we have that 

64f 2 

which is not monotone (Fig. |l|) nor, a fortiori, operator monotone. ("All stochastically monotone Riemannian metrics 
are characterized by means of operator monotone functions and . . . there exists a maximal and minimal among them" 
&) 
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l' ' ' '2' ' ' '3' ' ' '4' ' '5 t 
FIG. 1. Nonmonotone nature of the indicator function (^), corresponding to the probability distribution (^|) with u = —2 

D. The subcases m — 4, 5, 6 



Continuing on to the scenario m = 4, however, the legitimacy of this particular argument (sec. HC), concerning 
monotonicity, is undermined by the fact that the 16 x 16 analog of (^) and @ is not the same as that achieved using 
the probability distribution (^) with u — — 2 (or any other specific value of u). Nevertheless, these two 16 x 16 density 
matrices share the same zero-nonzero pattern and possess two similar sets of three distinct eigenvalues (^, ^ for 
the mean matrix based on (|^) setting u = — 2, and j^, Jg, for the one relying upon the measure of Zyczkowski 
et al) of multiplicities five, nine and two, respectively, corresponding to identical eigenspaces. (C. Krattenthaler has 
a formal demonstration that averaging over the Bloch sphere with respect to any spherically- symmetric probability 
distribution yields the same collection of eigenspaces for any to.) 

For the cases to = 5 and 6, we have found results of a similar nature. Using the U{2) measure of Zyczkowski et al, 
the eigenvalues for to — 5 are and ^ (with multiplicities six, sixteen and ten, respectively) and for to = 6, 

they are 57535, 2Tuo' TIBo anc ^ Tio ( W1 th multiplicities seven, twenty- five, twenty-seven and five, respectively). (These 
multiplicities are in conformity with the formula ( |To| ) below, reported in [po[ .) The corresponding averaged matrices 
and those based on ([}]) share the same zero-nonzero patterns, but do not fully match for any particular value of u. 
(We have not, however, explicitly confirmed that the eigenspaces are identical, as we anticipate.) 

It would be quite interesting, it would appear, for its possible utility in universal quantum coding ]2C| , |23[ |, to obtain 
a formula for general to for the entries and eigenvalues of the 2™ x 2 m averaged density matrix based on the measure 
of Zyczkowski et al. (Presumably, the eigenvectors are, for all to, those already given by Krattenthaler and Slater 
p0|.) Then, one could try to establish the asymptotics of the relative entropy with respect to this averaged matrix of 
the m-fold tensor products of the 2x2 density matrices. 



E. Explicit spin states 

For general to, in the framework of [^o| based on (||), the subspace spanned by all those eigenvectors associated 
with the same (<i-th) eigenvalue |2^, eq. (2.12)], 



1 r(|-u)r(2 + m-d-M)r(l + d-u) m 

2^ r(| + f - «)r(i - u) ' d-o,h---,l 2 



o,i,..., bd (9) 



corresponds to those explicit spin states sec. 7.5.j] |f-|2f] with d spins "up" or "down" (and the other to — al 
spins, of course, the reverse). (Eigenvectors — linearly independent, but not orthonormalized — are enumerated by 
"ballot paths" and given by formula (2.14) of ^0|.) The 2 m -dimensional Hilbert space can be decomposed into the 
direct sum of carrier spaces of irreducible representations of SU(2) x S m . The multiplicities of the eigenvalues |2(], 
eq. (2.13)], 
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fm - 2d + l) 2 (m + 1\ ,m, , N 

M ^= (m+1) ( d )' rf - ' 1 '-L2J ( 10 ) 

are the dimensions of the corresponding irreps. Thus, for m = 4, tt = —2, using and (|Io|), we have the results 
previously mentioned, A 40 = ^,A 41 = ^,A 4 2 = ^j, and M 40 = 5,Af 41 = 9, Af 4 2 = 2. Employing these two 
formulas again, we can easily ascertain that the eigenvalues of the 4x4 matrix (||) are (threefold) and | (unrepeated). 
For the 8x8 matrix (^]), we have the fourfold eigenvalue | associated with the orthonormalized eigenvectors, 

(0,0,0,0,0,0,0,1), (0,0,0,-^,0,-^,-^,0), ( ,-L,-L,0,i, 0,0,0), (1,0,0,0,0,0,0,0) (11) 

and the fourfold eigenvalue yj, corresponding to 

(0,0,0,--L,0,0,i,0), ( O ,0,0,-i,0,^,-i,0), (12) 

(o,--L,o,o,-L,o,o,o), ( o,--L,|,o,--L,o,o,o). 

Our investigation will now move on to cases for which N > 2. However, in this regard, let us point out that in 
their discussion of explicit spin-j states, Biedenharn and Louck |l7], p. 423] have remarked that "Unfortunately, for 
arbitrary n and k, the construction of the basis vectors . . . has never been given (fully explicitly), principally because 
of unsolved problems relating to the additional labels (a), which are required to specify a basis of the space v£ n \ 
and which imply a multiplicity of occurrence of the irrep ... of SU(2) x S n in the representation. However, for the 
case k = |, it may be proved, using standard character formulas . . . , that the representation . . . contains no multiply 
occurring irreps of SU(2) x S n ... For the case of the coupling of n spin-^ angular momenta, the indices (a) are 
not required in the notation ... for the basis vectors." ("The so-called missing label problem is a recurring one in the 
theory of group and Lie algebra representation theory. To be specific, suppose that A labels an irrep of a group G 
and v indexes a basis for an irrep k of a subgroup H C G. If the k irreps of H that occur in the space of the A irreps 
of G arc multiplicity free, then a basis {|A«i/)} is uniquely defined by the (Xkv) labels. However, if some k irrep has 
a multiple occurrence, an additional multiplicity index a is needed to label a complete basis {|AaKi/)}. This is the 
generic situation" || (cf. Q).) 

In a related context, Katriel, Paldus and Pauncz |3l) in an article entitled, " Generalized Dirac Identities and 
Explicit Relations between the Permutational Symmetry and the Spin Operators for Systems of Identical Particles," 
wrote that, "the most interesting problem to be considered is associated with the breakdown of the one-to-one 
correspondence between the total spin and the irreducible representations of the symmetric group for systems of 
identical particles with an elementary spin a > \ . . . For two particles the symmetric group is Si with only two 
irreducible representations, [2] and [l 2 ]. For a = i these two representations fully characterize the two possible total 
spin states 5 = and S = 1, respectively. However, for two a = 1 particles three different total spin states are possible, 
two of which are symmetric ([2]) and one of which is antisymmetric ([l 2 ])- This is the first example where we no longer 
have a one-to-one correspondence between the total spin states and the irreducible representations of the symmetric 
group. For three a = 1 particles the total spin is already insufficient to specify the irreducible representation; there 
are two different S = 1 states, each belonging to a different irreducible representation. This situation is repeated, 
more frequently, for states of four and five particles, but it is only for six a — 1 particles that an even more severe 
labeling problem appears: namely, two different sets of functions with the same total spin (S = 2) correspond to the 
same irreducible representation, [4,2]. For a — | the same labeling difficulties are encountered. Already for four 
such particles, more than one state with a given total spin corresponds to the same irreducible representation. For 
cr = 2 this degeneracy is present for three particles already. In other words, the representation of the symmetric group 
generated by the spin eigenfunctions is a reducible one, and the irreducible components can have multiplicities larger 
than 1" (cf. H-EE3). 

The abstract to this paper of Katriel, Paldus and Pauncz [|l] reads: "The well-known one-to-one correspondence 
between the eigenstates of the total spin for a system of spin-^ particles and irreducible representations of the 
symmetric group with up to two rows in the Young shape is the basis of interesting formal developments in quantum 
chemistry and in the theory of magnetism. As an explicit manifestation of this correspondence the class operators of the 
symmetric group are demonstrated to be expressible in terms of the total spin operator. This correspondence does not 
hold for higher elementary spins. The extension to arbitrary spin is investigated using Schrodinger's generalization 
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of the Dirac identity, which expresses the transposition in terms of two-particle spin operators. It is shown that 
additional operators, which for <J = \ reduce to the total spin operator, are needed for a complete classification. Some 
aspects of the formalism are developed in detail for a = 1. In this case a classification identical with that provided 
by the irreducible representations of the symmetric group is obtained in terms of the eigenstates of two commuting 
operators, one of which is the total spin operator." 

It appears that the research reported below is somewhat similar in nature to that of Katriel, Paldus, and Pauncz, 
but with the focus now on the eigenstates not of the N x N operators O k (in their notation), where eq. (38)] 
(s, = (s x i, s y i, s z i) being a one-particle spin operator) 

N 

O k = O k (N) = J2( s ^ k > ( 13 ) 

i<j 

but on the eigenstates of the N m x N m mean density matrix obtained by averaging with respect to the (normalized) 
[/(iV)-invariant measure. ("We have not been able to work out the general commutation relation between O k and 
Oi for arbitrary k and I larger than 1, but we would like to conjecture that O n ;n — 1,2, ...,2a form a set of 
commuting operators for arbitrary a. The missing proof of this conjecture is the most obvious loose end in the 
present investigation" |f3l|.) 




III. THE CASE TV = 3 
A. The subcase m — 1 

We have performed a series of computations in the framework of Zyczkowski et al Q for the case N = 3, utilizing 
the Euler-angle parameterization of S*C/(3), along with the invariant volume element, recently given by Byrd [ p5[ eq. 
(1) and p. 14] (cf. |3(], eq. (17)] and |37|, p. 5]). Averaging the 3x3 density matrices accordingly, we obtained the 
simple intuitive result (cf. (g)), 

(\ °\ 

(14) 

We have not, for the 3x3 density matrices been able to perform a strictly comparable averaging, utilizing the 
normalized volume elements of the maximal and monotone metrics, but certain interesting results based on them 
have, nevertheless, been obtained |6[|^]. (We first note, analogously to the 2x2 density matrices, that the integral 
over the entire eight-dimensional convex set of the volume element of the maximal monotone metric diverges. For the 
minimal monotone metric, on the other hand, such divergence appears not to occur. However, it is more problematical 
in nature, in that it appears not possible to symbolically integrate the volume element completely over the eight 
dimensions In particular, in ||, using a doMWe-limiting argument (necessitated by the non-normalizability of 

the volume element of the maximal monotone metric), a sii-dimensional marginal probability distribution was found 
based on a set of eight (separable) variables parameterizing the 3x3 density matrices. Its two-dimensional marginal 
probability over the simplex spanned by the diagonal entries (a, 6, c) was 

15(1 ~ a)VE . (15) 
4irVby/c 

The original symmetry between these three variables was broken by a series of transformations — suggested by work 
of Bloore J|^] — employed to obtain the eight variables, the separation of which was required in order to perform 
the necessary integrations to obtain (|l5|). It is, then, natural to associate the variable a with the middle of the three 
levels (the one inaccessible to a spin-1 photon, due to its masslessness) . The expected values with respect to ([l5]) are 
(a) = I, (6) = (c) = |. We, of course, note that neither of these values equals ^ as in jl^), based on the approach of 
Zyczkowski et al [Q. 
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B. The subcase m = 2 



We have also conducted the same form of averaging as used to get ( |14| ) for the 9x9 density matrices (the "two-trit" 
case |3l] [II]]) obtained by taking the twofold tensor products (m = 2) of the 3x3 density matrices (N = 3). The 
result obtained is the symmetric matrix (the £/(3) analog of the 4x4 matrix 
423])) 
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(16) 



where 5 = w .000368414. (The off-diagonal entries, lacking the constant 7r, correspond to the expected value 
of the product of symmetrically located off-diagonal entries — a^a^ — in the underlying 3x3 density matrix [oy].) 
Tracing the 9x9 density matrix (Rq) over one of the two constituent subsystems, we obtain the 3x3 diagonal density 
matrix (jl4|) with entries equal to IT 

Of the nine eigenvalues of (16), three are two are |, and the other four (unrepeated) ones can be paired as 



A fi = - 



7 



1296tt\/2 



.126216, 



^7 = 77 



7 



1296tt\/2 



.123784, 



and 



(17) 



A, = i 



/331 



.12816, 



An 



/331 



8 ' 12967T\/2 ' ' 8 12967rV2 

The three-dimensional subspace is spanned by the orthonormal set of eigenvectors 

1 



(o,o, 0,0,0, --L 0,-^,0), 



(0,0, --^=,0,0, 0,-^,0,0), 



(0, 



V2 



.12184. 



0,-^=, 0,0, 0,0,0), 



while the two-dimensional subspace is generated by 

1 - 1 ,0,0,0,0, ' 9 



(0, 
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3^' ' 3a/2 
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= ,0, 
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7 331 3V331 3V331 V331 



= , 0,0,0, 
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3V331 



)• 



The remaining four one-dimensional eigenspaces (corresponding to the sequence of eigenvalues (|1 
viewed as providing "labels" to the eigenspaces) are given by 



(-i,o, 



\.o.i 



1 



1 



2V2' ' 2' 2V2' 2y/2' 2V2 



,0), 
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r,0,-- 



1 1 
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,0), 



(18) 
(19) 

which can be 
(20) 



13 



2V2' 



2^2' 2\/2' 2^2' 



=)• 



/33T' 2^2' 2^2' 2V2' > 

We note that the constant 7r, appearing throughout the 9x9 averaged matrix (|l6j) and in the four eigenvalues (p^), 
is absent — similarly to the parameter u in the 2 m x 2 m analyses of Krattenthaler and Slater [ 20) — from the set of 
nine orthonormalized eigenvectors (^), ( |Ts| ) and (^0|). Thus, we are able, by replacing tt in (^) by a free parameter 
(w), to obtain a one-parameter family of 9 x 9 matrices, all the members of which have the same sets of eigenvectors 



8 



(|l~8|), ( |l~9| ) and (|o|). (We might speculate that there exists a family of probability distributions, parameterized by v, 
which would give these 9x9 density matrices, as their expected values, as we now know is, in fact, the case for the 
particular value v = tt [ cf. sec. VII] |.) For \v\ > ^J^- ps .0794116, the 9x9 matrix has nonnegative eigenvalues, so it 

is, then, a density matrix. If one sets v to either |j^= ~ .238235 or g 7^ ps .091662, two of the unrepeated eigenvalues 

become and ^. As u tends to ±oo, the four unrepeated eigenvalues ( |l7| ) all approach A. So, in that limit, we 
obtain a sextet and an "antitriplet" jl2], p. 311]. The antitriplet is spanned by the vectors ([lj), while (replacing ( |l9| ) 
and (J20|) ) , the sextet is spanned by 



(0,0,0,0,0,0,0,0,1), 



(0,0,0,0,0-^,0,-^,0), 



(0,0,-^,0,0,0,-^,0,0), 



(21) 



(0,0,0,0,1,0,0,0,0), 



(0,i,0,i, 0,0, 0,0,0), 



(1,0,0,0,0,0,0,0,0). 



"We know that the product of two SU(3) triplets decomposes into a sextet and an antitriplet. The sextet is sym- 
metric, whereas the antitriplet is antisymmetric. Therefore, both representations can be contained with the baryon 
multiplet which has mixed symmetry, but only the sextet can be contained in the multiplet of baryon resonances 
which corresponds to a totally symmetric representation" |l2|, p. 388]. 

C. The subcase m = 3 

For the subcase m = 3 of N = 3, the diagonal entries of the 27 x 27 density matrix, obtained by averaging according 
to the measure of Zyczkowski et al B , are 



(a, (3, j3, 13, 13, 7, f3, 7, 0, 13, (3, 7, (3, a, (3, 7, /?, (3, (3, 7, (3, 7, (3, (3, (3, (3, a), 



(22) 



where a = ^ ps .0516667, (3 — fts .0366667,7 = ifjjg ~ -030833. We have been able to determine the ojff-diagonal 
entries (of which, 272 are zero), as well. For instance, the (1,3) entry along with twenty-three others have the value 
og T n . Of the non-zero off-diagonal entries, most were of such a form, that is -p-. However, there were also three 

86407T & ' ' (7T ' 

distinct rational numbers: (thirty-six occurrences),^^, (eighteen occurrences) and gg^ (twelve occurrences). If 
we denote the entries of the underlying 3x3 density matrix by a^- (i,j = 1,2,3), then the cells in the 27 x 27 
density matrix yielding ^jL corresponded to entries of the form auaijdji, those yielding corresponded to entries 
of the form aaajkakj, and those yielding were associated with cells of the type aijajkati (where i,j, k now refer 
to distinct values of 1, 2 or 3). 

Seven of the twenty-seven eigenvalues were expressible as rational numbers. These were 



— re .0166667 (isolated), 



7 31 

ps .0291667 (fourfold), 

240 \ h 600 



.0516667 (twofold). 



(23) 



Among the remaining eigenvalues, there were three pairs, corresponding to the roots (A) of the sixth-degree even 
polynomial, where y = (7 — 240A)tt, 



5504 - 317043y 2 
The eigenvalues (A) are explicitly expressible as 



4986360?/ 4 - 19131876y b = 0. 



(24) 



1 



240 720tt\/2 



.0288541, 



240 720rr\/2 



.0294793, 



(25) 



1 



240 40571-^2 



.0286109, — 



1 



240 405tt\/2 



.0297224, 



and 



43 



7 



240 6480ttV2 



.0289389, — 



43 



240 6480ttV2 



.0293994. 
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There were also eight isolated eigenvalues (A = .0495426, .0496514, .0500807, .0515902, .0517431, .0532526, .0536819, 
.0537907), corresponding to the roots of the eighth-degree even polynomial, where y = (31 — 600A)7r, 

383127080633857021 - 18544605647405907654y 2 + 4501753947892101744?/ 4 (26) 

-351843587054438400y 6 + 8926168066560000?/ 8 = 0. 
The smallest and largest of these eight isolated eigenvalues are 



31 _ V233718T ^ Q 31 ^2337181 ^ 

600 1 620007T V / 2 600 1620007r v / 2 

If we replace the constant tt by the variable v (as was done in the N = 3, m = 2 case of sec. Ill B| [cf. ( pTj ) ]), then 
the seven rational eigenvalues are unchanged and the other twenty can be obtained from the solutions (A) of ( J24| ) and 
( p6| ) using now y = (7 — 240A)u and y — (31 — 600A)u. In the limit v — > ±oo, the eigenvalues of the matrix degenerate 
to (isolated), (sixteenfold) and ^ (tenfold), in agreement with the cardinalities of the multiplets of SU(3) 
constructed from three particles jl2], exer. 9.5], |^3|, eq. (11.5.4)]. 

The eigenvector corresponding to the lowest-lying (isolated) eigenvalue (^j) is 



(0, 0, 0, 0, 0, -a, 0, a, 0, 0, 0, a, 0, 0, 0, -a, 0, 0, 0, -a, 0, a, 0, 0, 0, 0, 0), (28) 



where a = 



D. The subcase m 



We have ascertained that on the diagonal of the 81 x 81 density matrix averaged according to the measure of 
Zyczkowski et al [0, there are three occurrences of 3^ w .0233333, eighteen of ^ pa .0122222, twenty-four of 
ji^j pa .0141667 and thirty-six of pa .0102778. We have only so far, however, been able to determine a relatively 
small number of the off-diagonal entries. For example, the (l,3)-entry has the value 86 4Q 07r , the (2,3)-entry, — 54 QQ 07r , 
and the (ll,12)-entry, 



37 



108000TT ' 



IV. THE CASE N = 2 x 3 
A. The subcase m — 1 



Let us consider the tensor product of an arbitrary 2x2 density matrix and an arbitrary 3x3 density matrix (cf. 
PJ). Then, we average the 6x6 result (corresponding to the subcase m = 1) over the convex sets of two-dimensional 
and three-dimensional density matrices, using the U (2) x U (3) product measure, in the natural extension of the work 
of Zyczkowski et al Q . We obtain a diagonal matrix with entries g , which is simply the same as the tensor product 
of (^) and (pjj), evidence of the statistical independence of the constituent density matrices. 

B. The subcase m = 2 

Now, however, let us similarly average the twofold tensor products (m = 2) of these 6x6 density matrices with 
themselves. Then, the eigenvalues of the averaged 36 x 36 density matrix (966 of the 1296 entries being zero) are 
i fa .0138889 (multiplicity three), i ps .020833 (two), ^ pa .0231481 (nine), ^ pa .0347222 (six) and the remaining 
ones (all explicitly expressible as fractions involving square roots and tt) are .0203067 (one), .0206307 (one), .021036 
(one), .0213599 (one), .0338445 (three), .0343845 (three), .0350599 (three) and .0355999 (three). For instance, the last 
of these values is an approximation (cf. (|l7|)) to ^ + 23 5 3 ^^ - If, as previously, we replace the occurrences of tt in 
(the denominators of certain entries of) the averaged matrix by a parameter v and let v — > oo (or, equivalently, simply 
set the entries in question to zero), we obtain the eigenvalues ^ (multiplicity three), ^ (six), ^fg (nine) and 
(eighteen). For the first of the corresponding four eigenspaces, all the nonzero entries of the (three) orthonormalized 
spanning eigenvectors were ±i, for the second and third eigenspaces, ±i or i-^j, and for the (dominant) fourth, the 

entries were 1 or \ or -h=. 
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V. THE CASES JV = 2x3x2 AND N = 3 x 2 x 2 



Let us continue further along the lines of the immediately preceding analysis (sec. IV). We construct a 12 x 12 
density matrix by taking the ordered tensor product of arbitrary 2 x 2, 3 x 3 and 2x2 ones. Then, by taking the 
tensor product of the result with itself, we obtain a 144 x 144 density matrix, which we average with respect to the 
product measure for U(2) x U(3) x U(2). The mean density matrix obtained had 16,886 of its 20,736 entries, zero. 
Let us, first, report the structure of the eigenvalues of this matrix, if we replace the occurrences in its off-diagonal 
cells of 7T by a parameter v, which we then let go to infinity, as we have previously in this series of analyses. There 
are, then, six eigenspaces, corresponding to the eigenvalues (mulitplicity six), ^ (twenty-seven), — jtj (twelve) , Ygo 
(fifty-four), ^ (fifteen), i (thirty). 

Leaving the constant ir unaltered, we obtain eighteen, rather than six, eigenspaces. The multiplicities of the 
eigenvalues ^ and are the same as in the v — > oo analysis, while that of is reduced from twelve to four, 
from fifty-four to eighteen and A from thirty to ten. The manner in which these reductions in dimensionality take 
place can be seen by examining the orthogonal 64-dimensional Hilbert space. It is composed of four eigenspaces of 
dimension nine (corresponding to the eigenvalues i± — 7= and T5n± — V33i \ {out Q f dimension five (associated 

y i o to 180 29160tt\/2 180 29160tt\/2 ; ' v 

with the eigenvalues ^ ± 12 ^^^ an< ^ YE ^~ mschiyl) an< ^ ^ our °^ dimension two (corresponding to the eigenvalues 
^± 7 —r= and o77T±— *^=). Of course, it is easily noted that 4x 9 = 54- 18, 4 x 5 = 30- 10 and 4x 2 =12-4. 

240 38880tt\/2 240 38880tt\/2 ; j 

We have also conducted a parallel series of analyses for TV = 3 x 2 x 2, thus, ordering the three density matrices 
in the initial tensor product differently. We, of course, again obtain a 144 x 144 mean density matrix (having, once 
again, 16,886 of its 20,736 entries, equal to zero). In the limit, v — > oo, an eigenanalysis yielded precisely the same 
results as the corresponding analysis in the 2x3x2 case. And in fact, leaving the constant tt unaltered, we obtained 
the same set of eighteen distinct eigenvalues and associated multiplicities as in that case too. 



VI. THE CASE N = 4 



For this analysis, we rely upon the (Hurwitz/Euler- angle) parameterization, together with the accompanying Haar 
measure, for U(N), for general N, given by Zyczkowski and Kus jl5], eqs. (3.1) -(3.5)]. For the subcase m = 1, the 
averaged 4x4 diagonal matrix has non-zero entries equal to \. For the subcase m = 2, the averaged matrix had the 
diagonal entries 

(a, /?, 7, k, (3, a, 7, n, 7, 7, e, C, re, k, £ 77), (29) 
where a = 2327 8 = 3947 -7 = 1759 K = 971 e = 1583 f = 2357 and n = 299 The onlv non-zero off-diaeonal 

wiieie ix 32400"" 64800'' 28800 ' ^ 17280 ' t 21600 ' H> 43200 dliu '' 3600' umjf uuu «au uu Uldguiliu 

entries (in the symmetric mean 16 x 16 density matrix) were: 6 7 g 7 Q in the (2,5)-cell; 2 3 g 9 Q in the (3,9) and (7,10)-cells; 
■^280 m the (4,13) and (8,14)-cells; and 4 j, 97 in the (12,15)-cell. We note the absence of the constant ir, in contrast 
to the scenarios of sees. QI, IV and 0, in which 3x3 density matrices were incorporated. 

There were seven distinct eigenvalues : 4r = .05 (having multiplicity six), 2 12 6 77 sa .0591204 (isolated), w 
.0623843 (two^_p^ w .071821 (three), w .0721528 (two), » .073287 (isolated) and ^ » .0830556 

(isolated) (cf. W2l cxer. 11.6, eq. (1)]). (The multiplets of a two-particle systems in the group SU(4) are a sextet 
and a decuplct [|42|, ex. 11.8(1)].) Orthonormal bases for the seven eigenspaces are easily constructed, in which each 
eigenvector has fourteen or fifteen components zero, and the others either equal to 1 or ±-^|. 



VII. THE CASE N = 5 



For the case, N — 5, m — 2, we have only so far been able to determine that the first two diagonal entries of the 
mean 25 x 25 density matrix are and y|g (which must also be the value of the sixth diagonal entry). We did in 
fact pursue a parallel (and more computationally manageable) analysis, in which the 5times5 unitary matrices were 
replaced by the 5x5 orthogonal matrices, and the appropriate Haar measure [[45"| , App. A], then, used. However, the 
mean density matrix for the case m = 1 was found to be null. 
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VIII. PARAMETERIZED FAMILIES OF N M x N MEAN DENSITY MATRICES 

In our previous analyses, we have employed the measures introduced by Zyczkowski et al [0 to find, for a specific 
choice of the dimension N and power m, a particular mean density matrix. However, it should be noted that the 
imposition by Zyczkowski et al of a uniform measure on the (N — l)-dimensional simplex spanned by the eigenvalues 
(ei,e2, . . • ,eAr) is somewhat arbitrary in nature. In fact, Bayesian principles would suggest that one might employ 
instead a Dirichlet distribution, 



T(N - q 1 - q 2 - . . . - q N ) _ gi _ 92 _ g 
T(l- qi )T(l- q2 )...T(l-q N ) ei 62 - e " 



(30) 



with its N parameters (qi,q 2 , ■ ■ ■ , <7jv) all set to \ jig, eq. (3.7)]. (The uniform distribution is obtained by setting 
the parameters all to zero.) This observation has led us to expand our analysis to include (symmetric) Dirichlet 
distributions with all N parameters set equal to q. (We also note that in the analysis of Krattenthaler and Slater |2(J 
a one-parameter (u) family of probability distributions (||) was utilized to obtain a one-parameter family of 2 m x 2 m 
mean density matrices.) 

A. N = 3,m = 2 



We have been able to implement this approach for the case N = 3, to = 2 of sec. [II B . All the members of the 
one-parameter family of 9 x 9 density matrices obtained possess the same zero-nonzero pattern as (|l^) . The diagonal 
entries | and ^ (corresponding, as noted, to q = 0) are replaced, now, by the more general entries, gpEqo an d 

i2^4"-3q) , respectively. The constant g is replaced by 2 i67r(4-3q) ano ~ *he off-diagonal entry, by 12 ( 4 "!_3 g ) . Instead of 
the three eigenvalues > we have ^ — g( 4 l 3g ) and in place of the double eigenvalue g , there is ^ + 18 ^_ 3q -^ ■ The two 
pairs of unrepeated eigenvalues ( |l7j) take the more general form, 



3-2 g 7 3-2 g 7 

6 6(4 - 3«) 324(4 - 3g)7rv^' ? 6(4 -3g) 324(4 - 3?)^' 1 J 

3-2q V331 . 3-2q 



6(4 -3g) 324(4 - 3g)7r\/2' 6(4 - 3g) 324(4- 3g)7iV2" 

The corresponding eigenspaces are precisely the same as those given by (|l8|)-(p0|). 

B. N = 4, to = 2 

We have pursued an analogous strategy for the N — 4, m = 2 analysis of sec. |v|, setting the four parameters of the 
Dirichlet distribution (p0|), again, all to q. The zero- nonzero pattern was the same as previously reported (that is, for 
the particular instance, q — 0). Then, the diagonal entries — as denoted in (|29| ) — took the form, 

_ 2327- 1620q _ 3947 - 3240g _ 1759 - 1440g _ 971 - 864g 

6480(5 - Aq) ' ^~ 12960(5 - Aq)' 7 ~~ 5760(5 - Aq) ' K ~ 3456(5 - Aq) ' ( ' 

_ 1583 - 1080g c _ 2357 - 2160? _ 299 - 180 9 

6 ~ 4320(5 - Aq) ' ^ ~ 8640(5 - Aq) ' 77 ~ 720(5 - Aq) ' 

The (2,5)-cell is now 12 mo(l-4q) , while the (3,9) and (7,10)-entries are ^f^^z^j- Also, the (4,13) and (8,14)-cells are 
3456 1 (°5-4 g ) and the ( 12 T5)-entry is 8640 1 ( 9 5 7 _ 4g) ■ 

The sixfold eigenvalue now takes the more general form, yg — 1& ^_ Aq ^ , the isolated eigenvalue 2 1 ieoo b ecomes 
JE ~ 4320(5-4q) ' wrme the twofold eigenvalue gUjj is transformed to yg — 1728(5-49) ' Additionally, the threefold 
eigenvalue 3 2 2 3 4 2 7 becomes ^ + 3240 1 ^ 1 _ 4(? ^ , the twofold eigenvalue changes to + 2 88o(5-4q) , the isolated eigenvalue 
2 1 1 5 g 8 o is converted to yg + 432o 2 (5-4q) ano - *he isolated eigenvalue g 2 ^ is transformed to yg + ggo^ig) ■ The associated 
eigenspaces were as reported above. 
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C. N = 4,m = 3 



We have begun to investigate this particular scenario and have established, among other items, that the (52,52)- 
entry of the 64 x 64 mean density matrix is 



the (52,61)-cell is 



while the (53,53)-entry is 



(33) 



(34) 



(35) 



43581 + 10<j(2160<7 - 6283) 
172800(5 -Aq)(3-2q) 

5026 - 2675 
172800(5 -Aq) (3 -2q)' 

58387 + 6g(6480g- 15979) 
311040(5 -4q)(3-2q) ' 

D. N = 2 

We have also for the case, N = 2, been able to obtain a iwo-parameter family of 2 m x 2 m mean density matrices 
by replacing the uniform measure of Zyczkowski et al by an (asymmetric) Dirichlet (or beta) distribution ([30]), with 
its two parameters q\ and qi not necessarily being equal. Nevertheless, the decompositions of 2 m -dimensional Hilbert 
space i nto t he associated eigenspaces of the mean density matrices appear to be of the very same nature as reported 
in sec. HE . If we do equate the two parameters (qi — q^ = q), then for the case, m — 2, we obtain the isolated 
eigenvalue 2 (3-2q) ano - ^ ne threefold eigenvalue gfen|^ ■ For to — 3, we have two quadruplets, gpEfg ancl Sppfg • 

Analogously to what was done in ]2lj for the 2™ x 2™ density matrices, one might investigate the nature of the 
Bures distance between two mean density matrices (for N > 2) corresponding to nonidentical sets of values of the N 
parameter(s) of the Dirichlet distribution (|30|). 

E. N > 4 

It seems quite natural to conjecture, based on the results of this section, that, in general, factors of the form 
(N + 1 — Nq) will appear in parallel results for N > 4. We remark, in this regard, that for the case N = 2, m = 4, the 
denominators of the three distinct eigenvalues of the 16 x 16 mean density matrix are all proportional to (5 — 2q)(3 — 2q) 
(cf. © - H)). 

IX. CONCLUDING REMARKS 

Of course, it would be of interest to obtain analogs of the results presented above for additional values of N 
and/or m, and to explore areas of further possible application (cf. |p6| , p7| ). The relation of the orthonormal bases 
of 7V m -dimensional Hilbert space that we have reported, to other bases, such as the Bell (N — 2, to = 2) and 
Greenberger-Horne-Zeilinger (N — 2, m > 2) ones, used in the processing of quantum information E^-H, merits 
investigation. (These latter bases correspond to maximally entangled states.) In this regard, we make the simple 
observation that in the N = 2,m = 2 case, the eigenvectors found by Krattenthaler and Slater |2Q] separate into a 
triplet (A/2,0 = 3) and singlet (M^i = 1), as with the Bell basis. The triplet (and, of course, the singlet) is known 
to be irreducible |42, p. 103]. We have, however, not formally established the irreducibility of the eigenspaces for the 
N > 3 cases reported above (sees. HI, ^ and VI). 



It would be desirable to better understand the significance of the process (selection rule) by which we have an- 
nihilated those off-diagonal entries in the averaged density matrices (based on tensor products incorporating 3x3 
density matrices) which are equal to a rational number divided by 7r (leaving undisturbed those off-diagonal entries 
which are simply rational numbers) and thereby obtained eigenvalues, the multiplicities of which correspond to the 
(permuationally-symmetrized) multiplets. In this regard, it would be of interest to determine the relations of the 
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approach taken here to the alternative one of Katriel, Paldus and Pauncz |nj (from whose article we have excerpted 



certain passages at the end of sec. [IE) 



We should also note for the benefit of those who might desire to further pursue the lines of investigation followed 
here that our results have, in large part, been based on MATHEMATICA |5(]] computations, in particular, multiple 
integrations over multivariate polynomials in trigonometric functions (sines and cosines of single angles). However, 
since MATHEMATICA would not directly integrate the cumbersome expressions (at least in an acceptable amount 
of time), we found it necessary at each integration stage to reduce the problem to its simplest possible (univariate) 
form, saving auxiliarly the unused information for the next integration step. In fact, it appeared necessary at several 
points to pursue even more subtle strategies in order to reduce the computational (space and time) burden. In any 
case, we have attempted to present the most extensive analyses within our capabilities of achieving. We should also 
observe that we have not investigated, to any major extent, the possible role for numerical, as opposed to exact or 
symbolic, integration methods. 
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